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Abstract 

In six male anesthetized, tracheotomized, and mechanically ventilated rabbits, 
we imaged subpleural alveoli under microscopic view (60 x) through a "pleural 
window" obtained by stripping the endothoracic fascia and leaving the parietal 
pleura intact. Three different imaging scale levels were identified for the analy- 
sis on increasing stepwise local distending pressure (_P ld ) up to 16.5 cmH 2 0: 
alveoli, alveolar cluster, and whole image field. Alveolar profiles were manually 
traced, clusters of alveoli of similar size were identified through a contiguity- 
constrained hierarchical agglomerative clustering analysis and alveolar surface 
density (ASD) was estimated as the percentage of air on the whole image field. 
Alveolar area distributions were remarkably right-skewed and showed an 
increase in median value with a large topology-independent heterogeneity on 
increasing P ld . Modeling of alveolar area distributions on increasing P ld led to 
hypothesize that absolute alveolar compliance (change in surface area over 
change in P^) increases fairly linearly with increasing initial alveolar size, the 
corollary of this assumption being a constant specific compliance. Clusters 
were reciprocally interweaved due to their highly variable complex shapes. ASD 
was found to increase with a small coefficient of variation (CV <25%) with 
increasing P ld . The CV of lung volume at each transpulmonary pressure was 
further decreased (about 6%). The results of the study suggest that the consid- 
erable heterogeneity of alveolar size and of the corresponding alveolar mechan- 
ical behavior are homogenously distributed, resulting in a substantially 
homogenous mechanical behavior of lung units and whole organ. 



Introduction 

Lung parenchyma is considered as a system composed by 
polyhedral alveoli opening on the lumen of alveolar ducts 
(Wilson and Bachofen 1982). The tissue structure sur- 
rounding alveoli is a complex interconnected network of 
macromolecules with viscoelastic and elastic properties. 
Alveoli also present a liquid interface with air implying 
the existence of surface forces. Essentially, whole mechani- 
cal lung behavior reflects the mechanical properties of 
single alveoli whose distending pressure is given by the 
sum of an elastic component, developed by the parenchy- 
mal structures, and a pressure component reflecting 
surface forces. Several studies have considered how 



considerable perturbation of physiological condition at 
lung periphery affect the intrinsic viscoelastic properties 
and/or surface forces creating regional heterogeneities in 
the mechanical behavior. Conditions studied include 
remarkable decrease in airways caliper following broncho- 
constriction (Tgavalekos et al. 2005), chemical destruction 
of the elastic component (Suki et al. 2003), and acute 
lung injury (Allen et al. 2005). A recent article by Wilson 
(2013) proposed that, in physiological conditions, regio- 
nal differences in volume oscillations may actually result 
from the complex interaction between the stiffness of the 
elastic elements of the alveoli and surface tension relaxa- 
tion. We wished to provide a contribution to this hypoth- 
esis, addressing specifically two related questions: (1) 
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assess differences in alveolar absolute and specific compli- 
ance: the former is defined as the change in volume over 
change in distending pressure, while the latter is defined 
as the absolute compliance over initial alveolar volume 
and represents the intrinsic mechanical property of the 
elastic components; (2) how the potential heterogeneity 
in alveolar mechanical properties impacts on whole lung 
mechanics. To this aim, we developed an experimental 
method that allowed us to image subpleural airway termi- 
nal units keeping pleural space intact; this approach 
allowed a clear and neat view of unrestrained subpleural 
structures exposed to changing physiological local dis- 
tending pressure. We estimated alveolar compliance under 
static conditions, thus neglecting a phase shift between 
elastic and surface forces, as hypothesized by Wilson 
(2013). 

Materials and Methods 

Animal preparation 

A general consensus for the experimental procedures used 
in our research activity was obtained from Milano-Bicoc- 
ca University Ethical Committee. Experiments were per- 
formed on six adult (New Zealand White) male rabbits 
(weight range 1—1.5 kg) anesthetized with a bolus of 
2.5 mL/kg of a saline solution containing 0.25 g/mL of 
urethane plus 10 mg/mL of pentobarbital sodium injected 
into an ear vein. The absence of an eyelid closure reflex 
was used as an estimate of the proper level of anesthesia 
during the experiment. When necessary, subsequent doses 
of 0.5 mL of anesthetic were administered through a 
venous line. Tracheostomy was performed, and animals 
were intubated. Before connecting the rabbit to the venti- 
lator, paralysis was accomplished by pancuronium bro- 
mide ( 1 mg/kg body weight initial dose, supplemented by 
0.33 mg/kg every 40 min). Mechanical ventilation was 
provided with a tidal volume (VT) of ~20 mL and a fre- 
quency of 16 breaths/min. The inspiratory flow signal was 
measured by a Fleisch-type pneumotachograph (Sensym, 
Milpitas, CA) and integrated in order to obtain lung vol- 
ume; volume drift was avoided by correcting flow offset. 
Tracheal and esophageal (through esophageal balloon) 
pressures were measured by pressure transducers. Lung 
volume values were normalized to their maximum value, 
assumed at 20 cmH 2 0. 

An "intact pleural window" was prepared in the sixth 
intercostal space to allow a view of the lower lobe where 
the cardiac artifacts are relatively small. The skin and 
superficial muscles on the right side of the chest were 
resected and, when reaching the layer of the intercostal 
muscles (indicated by "b" in Fig. 1A), a surface area of 
about 0.5 cm 2 was freed from muscles down to the 
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Figure 1. Experimental preparation. (A) Macrophotography on 
creating a "pleural window": (a) intact chest wall, (b) intercostal 
muscle layer, (c) denuded portion of parietal pleura ("pleural 
window"). (B) Part of a pleural window under the microscopic view 
with a magnification of 60x: alveoli are clearly visible where the 
endothoracic fascia was stripped off the denuded parietal pleura 
(c1), as compared to the unstripped portion (c2). (C) An example of 
manually segmented alveolar units (digital zoom of a 60x 
magnified microscopic image) at P u = 4.5 cmH 2 0. Four 
nonoverlapped ROI are shown together with the corresponding 
percentage values of alveolar surface (alveolar surface density, 
ASD). (D) The same alveolar population of (C) with manually traced 
alveolar clusters (black lines) obtained by the contiguity-constrained 
hierarchical clustering analysis. Alveoli not belonging to any clusters 
are not encircled. 

endothoracic fascia. Under stereomicroscopic view and 
with fine forceps, the endothoracic fascia was then carefully 
stripped, exposing the transparent parietal pleura (about 
10 fim thick) through which freely moving subpleural alve- 
olar structures could be neatly visualized (area indicated by 
"c" in Fig. 1A). The microscopic image in Figure IB high- 
lights the difference between the stripped area (indicated 
by "cl") and the unstripped one (indicated by "c2"). 

Experimental protocol 
Image acquisition 

The animal was placed in left lateral decubitus under the 
field of a fixed microscope (Nikon's Eclipse FN1, Nikon 
Corporation Instruments Company, Tokyo, Japan) con- 
nected to a video camera (CoolSNAP B2 , Photometries, 
Tucson, AZ) which was interfaced through a IEEE- 1394 
data-transfer interface card with a personal computer 
equipped with an image-processing software (MetaMorph® 
System, Molecular Devices, Sunnyvale, CA). Total 
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microscopic magnification was 60 x. A LED ring-light 
illuminator was anchored to microscopic optics to provide 
a uniform lighting of the alveolar field. The depth of field 
of the objective used (4x, NA = 0.1) was about 
50 microns. As we focused the visceral pleura, the micro- 
scopic image referred to a depth of 50 microns within the 
lung parenchyma. We also refer to common histology 
(hematoxilin-eosin) performed in our laboratory on sagit- 
tal sections of the lung fixed at a transpulmonary pressure 
corresponding to functional residual capacity (FRC) to esti- 
mate the average thickness of the visceral pleura (Rivolta 
et al. 2011). 

We imaged the alveolar units of the right lower lobe, 
where no interlobar fissures crossed the optical field. 
Alveolar units included both alveoli and alveolar sacs. The 
change in alveolar size was assessed during two/three 
recruitment-derecruitment maneuvers: mechanical venti- 
lation was stopped at the end of expiration and lung was 
recruited from FRC up to an alveolar pressure (-P a i v ) of 
20 cmH 2 0 (100% of inspiratory capacity, IC) by injecting 
air through a syringe and controlling intratracheal pressure. 
Each step involved an increase/decrease of 4 cmH 2 0. 
During the maneuver, an image was acquired at each pres- 
sure step, when quasi-static conditions were reached (about 
2 sec for each pressure level). The displacement of the 
alveolar field under microscopic view required a manual 
focus adjustment at each pressure step. At the end of the 
experiment, the animals were euthanized with an anesthetic 
overdose. 

Morphological and morphometric analysis of 
alveolar units 

For each image, alveolar units (both alveoli and alveolar 
sacs) were manually segmented along the sharp gradient 
in optical density between the gas phase and the intersti- 
tial tissue in order to derive individual alveolar surface 
area (Fig. 1C) and to obtain then their frequency distri- 
butions. The latters did not include alveolar sacs, whose 
geometry depart from circularity. Given a rectangular 
region of interest (ROI) on the image (with a variable size 
around 200 x 300 micron, Fig. 1C), the alveolar surface 
density (ASD) was estimated as the ratio of alveolar sur- 
face area to the area of the ROI; up to four equally sized, 
minimally overlapping, randomly selected ROIs contain- 
ing at least 10 alveoli were analyzed for each image to 
obtain an average value for ASD at each alveolar pressure. 

Local distending pressure 

Given the animal position, the imaged lung units of the 
right lower lobe were placed in the less-dependent por- 
tion of the pleural cavity, that is, at 100% of lung height. 



For these units, it was possible to obtain an estimate of 
the local distending pressure (P ld ) defined as: 

Pid = Pah — Pips 

where Pi ps is the local pleural pressure. P lps was derived from 
figure 7 of Agostoni and Miserocchi (1970), a study that 
reported direct measurements performed at different lung 
heights and for different P a i v values in the same species. 

Cluster and fractal analysis 

Cluster analysis was carried out to describe the shape and 
topological distribution of clusters of alveoli having simi- 
lar area values and was based on contiguity-constrained 
hierarchical agglomerative clustering (Recchia 2010). This 
analysis allows to group alveoli based on two inclusion 
principles: the first one is based on similarity in surface 
area, the second is contiguity. 

The contiguity constraint is initially computed on single 
alveoli that are considered adjacent if their centroids are 
not farthest than two times their average surface diameter 
and no other alveoli are found on the line joining their 
centroids. A representative example of this clustering 
analysis is presented in Figure ID, where alveolar clusters 
are manually encircled by black lines. Clusters were geo- 
metrically characterized by computing their circularity, 
defined as 4n x (area/perimeter 2 ) grouping data at _P 0 . 

Statistical analysis 

Kolmogorov-Smirnov normality test was performed on 
alveolar populations using Matlab®. Spatial statistical 
analysis (semivariogram) of areas of single alveoli was 
performed to express the relationship of the respective 
variances as a function of increasing distance ("lag") 
between individual alveoli, in order to assess whether 
there is a distance dependence for alveolar area distribu- 
tion. Logarithms of area values were used to obtain nor- 
mal distributions, on which semivariogram analysis was 
carried out, using geoR package (R, version 3.0.1). Semi- 
variogram could be computed only up to a maximum 
distance corresponding to half the window size (about 
500 Aim). 

Results 

Data will be presented in accord to three scale levels: 
alveoli, cluster of alveoli, and overall image/whole lung. 

Alveolar level 

Representative distributions of alveolar areas pooled from 
all rabbits are presented in Figure 2A and B for 
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Figure 2. (A and B) Distributions of alveolar areas at Pu = 3 and 16.5 cmH 2 0 (P 0 and respectively), obtained by pooling data from all 
rabbits. (C) Box plot describing the alveolar surface areas from all rabbits at different levels of alveolar pressure (only the inflation limb is 
shown). For each pressure, the gray box shows a 25-75% percentile range; the median value is represented by the horizontal bar. The 
whiskers above and below the gray boxes encompass 80% of the variability in area values. (D) Semivariogram of alveolar areas at P 0 . 



P alv = 0 cmH 2 0 and P alv = 20 cmH 2 0, corresponding to 
local distending pressure P ld of 3 and 16.5 cmH 2 0, 
respectively referred, from this point onwards, as P 0 and 
Pi. The total number of alveoli analyzed for each image 
ranged between 100 and 150, providing an overall num- 
ber of alveoli in a range of 1000 at each pressure. These 
distributions were found to be right-skewed, that is smal- 
ler alveoli are more numerous than the greater ones. The 
semivariogram presented in Figure 2C shows that the var- 
iance of the alveolar areas at P 0 is independent of topo- 
logical distribution, as it oscillates by about ±3% around 
the variance of the alveolar population, named "sill." 
Similar results were obtained for the same analysis carried 
out at increased distending pressure, with an average 
oscillation of ±7% around the "sill." 

Figure 2D shows that the median values of alveolar 
areas, pooling data from all rabbits, increase with P ld , 
with a remarkably increasing wide variability (expressed 
by the 25fh-75th percentiles, dashed areas). Kolmogo- 



rov-Smirnov normality tests failed for distributions at 
all distending pressures. Also log-normality test, per- 
formed by applying normality test on the logarithm of 
data, failed for each distribution. The coefficient of vari- 
ation (CV) of each distribution was computed by divid- 
ing the interquartile range (75th percentile value minus 
25th percentile value) by the median value and did not 
show any specific trend to increase with increasing dis- 
tending pressure: for local distending pressure of 3, 4.5, 
8.5, 10.5, 13, and 16.5 cmH 2 0 on the inflation limb, 
the CV values were 88%, 79%, 87%, 92%, 79%, 80%, 
respectively. 

Cluster level 

Figure ID allows to appreciate the complex matching of 
clusters due to their highly variable shape. Figure 3 shows 
that cluster shapes are far from circularity and quite dis- 
persed. 
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Whole image level 

The average values of ASD ± SD, obtained by pooling 
values from all animals, is reported in Figure 4A and 
shows an increase up to a P ld of 10 cmH 2 0, with CV at 
each distending pressure not exceeding 25%. Figure 4B 
shows the overall lung volume versus distending pressure 
curve and allows to appreciate that the average CV of vol- 
ume is further reduced to about 6%. 



Discussion 

Previous attempts of imaging of subpleural alveoli to 
some extent largely interfered with unrestrained alveolar 
movement (Liu et al. 2011). For example, in Schwennin- 
ger et al. (2008) and Tabuchi et al. (2008), the visceral 
pleura was stuck to a transparent window. By using a 
tomographic technique as optical coherence tomography 
(Meissner et al. 2010) or confocal microscopy (Namati 



et al. 2008), a three-dimensional (3D) analysis could be 
performed. In this case, however, less pressure steps could 
be done in in vivo experiments due to the time required 
by the 3D acquisitions. 

The experimental approach used in this study allowed 
to investigate the morphology of the subpleural alveoli by 
maintaining pleural space mechanics intact. A notable 
limitation was that it was not possible to follow the mor- 
phological change in the same units. On the other side, 
the method had the advantage to allow unrestrained 
movement of subpleural units, exposed to physiological 
distending pressure, by preserving the integrity of the 
respiratory system. This allowed to relate the morphology 
of the alveoli, assessed through a two-dimensional analy- 
sis, to the expected local distending pressure (Agostoni 
and Miserocchi 1970). This last point is of remarkable 
importance as it is not available by other techniques used 
to investigate the morphology of individual alveoli within 
the lung parenchyma. Another potential limitation is 
represented by the fact that only subpleural units were 
studied and the discussion concerning their mechanical 
behavior does not allow us to extrapolate to the rest of 
the parenchyma. In fact, the areas of the subpleural alve- 
oli are controlled by the lines of attachment of alveolar 
walls to the pleura and the displacements of these attach- 
ments relate to the displacement of the pleural membrane 
which has been shown to expand isotropically as lung 
volume increases (Lehr et al. 1985). However, recalling 
data from Gil et al. (1979), no differences in surface area 
expansion have been described between subpleural and 
core alveoli from fixed rabbit lungs. Finally, our imaging 
system did not allow a 3D reconstruction. Given the 
depth of field of our objective (50 microns), we essen- 
tially imaged the terminal airways that may include both 
alveoli as well as alveolar sacs, as confirmed by sagittal 
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Figure 5. (A) Absolute compliance estimated by comparing subpopulations drawn from alveolar area distributions at P 0 and P, (continuous 
line). (B) Specific compliance derived from the absolute compliance, as detailed in the text (continuous line). Dashed line in (A) and (B) 
correspond to the assumption of a constant specific compliance. (C) Modeling of the experimental alveolar area distribution at P\ (white dotted 
points), assuming either a variable (black dashed line) or fixed (gray dashed line) specific compliance. 



histological sections from our lab; as specified in the 
methods, we only considered surface area values of the 
alveoli. One may wonder whether the observed heteroge- 
neity could in part be due to the fact that the focal plane 
cuts the alveoli at different positions with respect to their 
centroids. This would lead to an error in the measured 
areas, representing an underestimate relative to the area 
corresponding to the plane of the alveolar centroid. From 
the histological sections, we also estimated an average 
thickness of the visceral pleura of 10 ± 2 microns, so that 
all alveolar walls are essentially at equal distance from the 
pleural surface. Therefore, the underestimate would only 
occur for larger alveoli when the centroid is much below 
the focal field. This bias, however, is presumably constant 
at all pressures, so that it would not affect relative 
changes in alveolar areas observed in increasing distend- 
ing pressure. 

By considering the histograms of Figure 2A, we will 
attempt an interpretation of the change in distribution of 
the alveolar caliper with increasing Py to derive indica- 
tions on the mechanical properties of the alveoli. In par- 
ticular, we propose to interpret our data in terms of 
alveolar distensibility. The latter property was estimated 
by the ratio of the increase in alveolar surface area for an 
increase in alveolar distending pressure from P 0 to Pi, 



namely C abs = (A 1 - A 0 )/{P 1 — P 0 ), which represents an 
average index of alveolar absolute compliance in the range 
of pressure considered. We also defined specific alveolar 
compliance as C sp = C a b s /A 0 . As we cannot follow the 
same alveoli on increasing Py, we attempted to derive 
both absolute and specific alveolar compliance from the 
frequency distribution of alveolar areas at P 0 and Pit as 
detailed in the Appendix Al. We implicitly assume that, 
within a frequency distribution at a given Py, alveoli of 
different size are actually all exposed to Py, for mechanical 
equilibrium. Figure 5A shows that the computed estimate 
of C a bs (continuous line) increases sevenfold for a range 
of increase in A 0 including 90% of alveolar population. 
Correspondingly, one can appreciate that specific compli- 
ance C sp (Fig. 5B, continuous line) remains essentially 
steady (0.12 -=- 0.15 cmH 2 0 _1 ) for the same range of alve- 
olar population. To simplify the model, we assumed a 
constant value for specific compliance (see Appendix Al), 
reported by the dashed lines in Figure 5B. Obviously, a 
constant value for C sp implies a linear increase in C a b s 
(Fig. 5 A, dashed line). To validate this assumption, we 
show in Figure 5C that a similar satisfactory modeling of 
the experimental distribution (white points) at ri was 
found considering either the variable or fixed specific 
compliance (dashed black and gray line, respectively). One 
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shall comment that the derived absolute compliance 
values obviously reflect the shape of the experimental 
area distribution curves and the modeling assumption of 
the correspondence of the z-th subpopulations at various 
Pi* 

Consideration relative to the difference in 
the mechanical behavior of lung units 

Results of this study reveal that the considerable heteroge- 
neity in alveolar size (Fig. 2C) does not show a topologi- 
cal dependence, as demonstrated by the semivariogram of 
Figure 2D, at least over a limited range of distance 
(500 ;im). Furthermore, the clusters also appear to be 
interweaved (Fig. ID), having a highly variable shape 
(Fig. 3). Therefore, we assume that alveoli are not inde- 
pendent entities but are exposed to mechanical interde- 
pendence regardless of the fact that alveoli are or are not 
subtended by the same airway. 

Obviously, we are far from extending our consider- 
ations to the whole lung functions where inhomogeneous 
behavior is known to be present concerning various 
aspects such as gas mixing, regional volume change, due 
to several interacting factors, such as gravity, airway resis- 
tance, etc. Interestingly, scaling up from single alveoli to 
randomly selected alveolar regions, a relatively homoge- 
nous mechanical behavior with minimal hysteresis of 
overall alveolar expansion is observed, as suggested by the 
low CV of ASD (Fig. 4A). This occurs despite the fact that 
differences in alveolar morphology imply differences in 
absolute compliance (Fig. 5). Scaling up to the level of the 
whole lung, a relative homogenous behavior is also 
observed given the low CV of lung volumes at various 
transpulmonary pressures. Thus, to summarize, we mea- 
sured the CV of volumes at alveolar, ROI and whole lung 
levels and found that heterogeneity is maximal at the alve- 
olar level. Note that for the whole lung expansion, various 
structures with different mechanical properties become 
involved (alveolar ducts, bronchioles, bronchi). In fact, on 
comparing Figure 4A and B, one can deduce that above 
10 cmH 2 0 a contribution to increase in lung volume 
must come from lung structures other than the alveoli 
under analysis and this also implies an increase in hystere- 
sis from ASD to whole organ, confirming previous find- 
ings obtained on isolated lung lobes using a synchronized 
stroboscopic photography (Lehr et al. 1985). In conclu- 
sion, to answer the questions addressed in the introduc- 
tion, our data suggest that the mechanical behavior found 
by scaling up to lung units and to the whole organ results 
from a non topological-dependent morphological and 
mechanical heterogeneity at microscale level (alveoli 
and cluster). The relative constancy of the estimated alve- 
olar specific compliance could be an important factor 
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to account for morphofunctional homogeneity on increas- 
ing scale level. Indeed, we might consider that, in an inter- 
dependent environment, a homogenous distribution of 
specific compliance of elastic components might reduce in- 
terregional differences in parenchymal forces during lung 
volume changes. Of course, this working hypothesis 
requires an experimental validation by following the 
expansion of the same alveoli on increasing distending 
pressure. Regarding the fact that no hysteresis was observed 
in alveolar areas, we can make the following comment: on 
one side, this may be due to the fact that we could actually 
not follow the same alveoli; on the other, we may suppose 
that hysteresis arises from other structures such as alveolar 
ducts. The hysteresis observed in lung volume change, 
despite a small CV of lung volume (Fig. 4B), may be 
ascribed to a multiscale effect of hysteretic factors on lung 
units having different mechanical properties. This can be 
the object of a further study requiring the use of imaging 
techniques allowing to explore internal parenchymal 
structures. 
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Appendix 

We defined subpopulations from the frequency distribu- 
tions at P 0 and Pi by dividing the overall populations in 
50-quantiles and we assumed that each i-th subpopula- 
tion at P 0 fell in the corresponding i-th subpopulation at 
Pi. The absolute compliance of the i-th subpopulation, 
Qbsf) from P 0 to Pi was estimated from the increase in 
median value of the subpopulation itself (Ai. — Ao f )> 
divided by the increase in P ld : 



P1-P0 

The specific compliance index was given by 



(i) 



(2) 



A constant value of specific compliance was estimated 



by 



An 



A 0 



(3) 



(A 0 where Ai and A 0 are the median value of distribu- 
tions at P 0 and Pj ) 
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